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We examine the gluon propagator in the Laplacian gauge in quenched lattice QCD as a function 
of the number of colours. We observe a weak dependence on Nc over the whole momentum range. 
This implies an almost A^c-independent gluon pole mass in units of the string tension. 
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I. INTRODUCTION 

Lattice QCD provides a well suited theoretical frame- 
work for the study of nonperturbative quantities as a 
function of the number of colours. Although experimen- 
tally inaccessible, QCD for — 2 and A^^ = 4 is in- 
teresting from the theoretical point of view. Recently 
SU{Nc) Yang-Mills theories were investigated for vari- 
ous Nc by two groups with emphasis on the eval- 
uation of string tensions and glueball masses and their 
large Nc limit. It is the purpose of the present work to 
examine the gluon propagator for different numbers of 
colours. Like in our previous work in SU{3), we evalu- 
ate the gluon propagator in the Laplacian gauge which, 
unlike the more usual Landau gauge, is free of spurious 
Gribov copies on the lattice. Therefore no doubt can be 
cast on our results due to the lattice gauge fixing pro- 
cedure. Furthermore theoretical arguments have been 
presented supporting gauge invariance of the pole mass 
of the transverse gluon propagator j|] and it is therefore 
an interesting quantity to investigate as a function of Nc 
even in a particular gauge. Studying the dependence of 
the gluon mass on Nc has theoretical implications for 
confinement and chiral symmetry breaking. It has been 
argued that a non-zero gluon mass is connected to vor- 
tex condensation and to the glueball mass Since the 
Nc dependence of the glueballs has been investigated, a 
study of the gluon mass can shed light on their relation. 
In Dyson Schwinger studies the form of the gluon propa- 
gator affects that of the quark propagator, which in turn 
provides a determination of the quark mass and therefore 
is connected to chiral symmetry breaking. 

In previous studies we considered the gluon prop- 
agator in SU{3), and presented evidence for a pole which 
survives the continuum limit. In this work we present re- 



sults in SU{2), SU{3) and S'?7(4) at similar couphngs. 
We refer the reader to rcf. js) for our notation and the 
details of our approach. 



II. LAPLACIAN GAUGE FIXING 

The Laplacian gauge consists of rotating along a fixed 
orientation the local color frame built from the lowest- 
lying eigenvectors of the covariant Laplacian. It is well 
defined in the continuum theory but non-perturbative. 
On the lattice, it has the virtue of being unambiguous 
(except for genuine Gribov copies arising from accidental 
degeneracy of Laplacian eigenvalues), unlike the lattice 
Landau gauge. 

In the case of SU{2), the Laplacian eigenvalues are 
twofold degenerate due to charge conjugation symme- 
try: if / is an eigenvector so is a2f*. At each lattice 
site, one can construct from these two vectors a 2 x 2 
matrix which is then projected onto SU (2) and defines 
the required gauge rotation. 

In SU{3) and SU{i), we identify the Nc-1 lowest-lying 
eigenvectors of the Laplacian. From these eigenvectors 
we construct at each lattice site an SU{Nc) matrix by 
using a Gram-Schmidt orthogonalization, as described in 
detail in ref. [||. 



III. RESULTS 

In order to make the comparison of our results in 
SU{2), SU{3) and SU{A) easier, we consider /3 val- 
ues which yield similar lattice spacings in units of the 
string tension a. We take /3 = 2.4, 6.0 and 10.9 for 
SU{2), SU{3) and SU{4) respectively, which give o^ct = 



1 



0.2634(16) §, 0.2265(55) § and 0.2429(14) (| for the 
three gauge groups. We can then express our SU{Nc) 
results in "physical" units using y/a = 440 MeV. The in- 
verse lattice spacing is 1.67(1), 1.94(5) and 1.81(1) GeV 
for SU{2), SU{3) and SU{4) respectively. 

We first examine SU{2). For this gauge group, several 
lattice studies of the gluon propagator in the Landau and 
Coloumb gauges exist |]l0| , |lT| . We present here the first 
study in the Laplacian gauge where no issues of lattice 
Gribov copies arise. It is therefore important to have lat- 
tice artifacts under control. The volume dependence is 
investigated by considering the SU{2) gluon propagator 
for lattices of size 4^, 8^, 16"*, 16^ x 64 and 32^. The 
quantity of interest is the transverse part, D{q^), of the 
gluon propagator. This is shown in Fig. 1 in lattice units 
for the various volumes after making a cylindrical cut | [T^ ] 
in the momenta which eliminates the leading lattice ar- 
tifacts due to lack of rotational symmetry. For non-zero 
momenta, we subtract the longitudinal part of the propa- 
gator to project out its transverse part. For zero momen- 
tum this subtraction is not well defined, and we keep the 
complete propagator. We include results from the very 
asymmetric 16^ x 64 lattice. Even in that case, only small 
finite-size effects are visible in the infrared (momenta less 
than ~ 0.5 GeV in "physical" units). 
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FIG. 1. The 5*17(2) gluon propagator in lattice units for 
various volumes. 



A similar study of the volume dependence is carried 
out for the SU(4) gluon propagator for lattices of size 
A'^, 8^, 8^ X 32 and 16^ x 32 and the results are shown 
in Fig. 2. As for SU(2) we observe only small finite-size 
effects in the infrared for the very asymmetric 8'^ x 32 
lattice. 
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FIG. 2. The SU{4:) gluon propagator in lattice units for 
various volumes. 



One conclusion to be drawn from Figs. 1 and 2 is that 
the zero momentum gluon propagator has a finite value 
as the volume is increased. In fact for lattices of spatial 
size greater than 0.8 fm it reaches a constant value. This 
behaviour also holds for SU{3) |^,^. 

In Fig. 3 we make a comparison of the renormalized 
zero momentum gluon propagator, I'(O) in "physical" 
units, collecting all our results for the three gauge groups. 
To relate the bare lattice propagator to the renormalized 
continuum propagator I?r,((7; /i) one needs the renormal- 
ization constant Z3(/i, a): 



a^D{qa) = Z^ip,, a)Dji{q; IJ.) 



We take as our renormalization condition 



(1) 



(2) 



at a renormalization scale fi which allows a de- 
termination of Z^[^,^a). Taking the renormaliza- 
tion point at /i = 1.94 GeV we find for SU(3) 
Z3^^^^^(^,a) = 2.51(1), and for SU(2) and SU(4) 
the ratios zl^^^\^i,a) / Z^^^'^\^i,a) = 1.07(4) and 

Z^^^^\^ji,a)/Z^^^^\li,a) = 0.96(2) respectively. 

All the results for various volumes fall rather well on 
a universal curve demonstrating independence on Nc- In 
particular both in SU(2) and SU(3) they reach a plateau 
at approximately similar volumes. Since the zero mo- 
mentum propagator gives a measure of the range over 
which the gauge fields are correlated, these results show 
that the gauge fields decorrelate in SU{Nc) at a distance 
of ~ 0.8 fm or about 1.7 times the confining correlation 
length. 
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FIG. 3. Volume dependence of the renormalised zero mo- 
mentum gluon propagator ©(0) for SU{2) (crosses), SU{3) 
(triangles) and S[/(4) (stars). The dashed line is the best fit 
to the form aexp{~V/Vo) + c yielding Vb ~ 0.1 fm**. 



In Fig. 4 we plot the results for the renormalized trans- 
verse gluon propagator for lattice sizes for which we ex- 
pect volume effects to be negligible, namely 16'^ x 32 for 
SU{3) and SU{4), and 32^ for SU{2). Remarkably, the 
overall shape of the gluon propagator shows no clear Nc 
dependence. 
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FIG. 4. The momentum dependence of the renormalised 
transverse gluon propagator for SU(2) (crosses), SU(3) (tri- 
angles) and SU{4:) (stars) . 



However, infrared information is suppressed when 
showing q^D{q^). To make any low-momentum differ- 
ence clearly visible, we plot in Fig. 5 the SU{2) and 
SU{4) propagators versus the SU{'i) propagator for var- 
ious momenta. If there would be no iVc-dependence 



the results should fall on the diagonal y = x. Small, 
but systematic differences are visible between SU (3) and 

SU{2) S't/(4) where the propagators are the largest, 

i.e. in the infrared, for "physical" momenta smaller than 
~ 1 GeV. The SU (3) propagator is slightly larger, which 
will translate into a slightly smaller gluon pole mass. 
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FIG. 5. The renormalised transverse gluon propagator for 
SU{2) (crosses) and SU{4:) (squares) versus that for SU{3). 



To investigate more closely the q dependence of the 
gluon propagator we perform various fits to phenomeno- 
logical Ansatze. We fit the transverse propagator to 
Cornwall's Ansatz which allows for a dynamically gener- 
ated gluon mass and is physically motivated. We also fit 
to two parametrizations which appeared recently in the 
literature, namely Model A and the mass parametri- 
sation of ref . lliol : 
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1 + qym^ [log{l + q^/ml)]^^^'^\ 



(3) 



In Fig. 6 we show the quality of the fits to the three 
Ansatze: Cornwall's Ansatz has three fit parameters, 
Model A four and the mass Ansatz five. In accord with 
ref. which studied the SU{2) Landau gauge, we find 
that the mass Ansatz works well for SU{Nc). This is not 
so surprising since it has the largest number of parame- 
ters. We find no compelling reason to prefer this Ansatz 
to the other two. Cornwall's Ansatz is the only one to 
allow a smooth analytic continuation to —q^ and we will 
use it for the determination of the pole mass. 



3 




SU(3) 

Cornwall 
Model A 

mass 




12 3 4 5 

q (Gev; 

FIG. 6. Fits of the gluon propagator to Cornwall's Ansatz 
(solid line) , Model A (dashed line) and the mass parametri- 
sation (dash-dotted line), in SU{2) (upper), SU{3) (middle) 
and SU^i) (lower). 



The salient features of the infrared behaviour of the 
gluon propagator in the three gauge groups show up by 
plotting its inverse as a function of q^. From Fig. 6 it 
can be seen that the overall behaviour is very similar for 
the three gauge groups although, as shown in Fig. 5, not 
identical. To make this observation more quantitative we 
investigate the gluon pole mass by analytic continuation 
to negative values of as shown in Fig. 7. In addition 
to Cornwall's Ansatz we use polynomials of increasing 
degree in to test the model dependence of the analytic 
continuation. 
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and analytic 
^ for the determination of a pole mass: Up- 



FIG. 7 
continuation to —q 
per graph for SIJ{2), middle for SC/(3) and lower for SIJ{^). 
We show Cornwall's Ansatz (solid line), a linear (dotted line), 
a quadratic (dashed-dotted line) and a cubic (dashed line) 
polynomial in g^. 



Taking as a mean the value extracted from Cornwall's 
Ansatz and as a crude indication of the systematic error 
the results of the polynomial fit, we find for the pole mass 
(1.85 ± 0.30)Vct, (1.52 ± 0.06)v^ and (1.74 ± 0.27)v^ 
in S'C/(2), Sl}{i) and 5[/(4) respectively A study of 
the gluon propagator for SU{2) in Landau gauge |ic| ] 
found that the low momentum behaviour supported a 
non-zero mass of ~ (1.48±0.05)-y/o'. Since only the mass 
parametrisation was used in that study and no analytic 
continuation was performed, this value can be regarded 
as a rough estimate consistent with our pole mass. 
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IV. CONCLUSIONS 



We have extended a previous lattice study of the gluon 
propagator in the Laplacian gauge to SU (Nc). In SU{2), 
we have made an analysis of the finite-size effects, which 
indicates that the Laplacian gauge gluon propagator at 
zero momentum approaches a non-zero constant value in 
the infinite volume limit, as in SU{3). In fact, this finite 
non-zero value is independent of TVc within our statistics. 
We find V{0) = 5.85(5) GeV-^. 

The overall dependence of the gluon propagator on Nc 
appears extremely weak, and is barely visible, even for 
infrared momenta, even for Nc = 2. For all Nc = 2,3,4, 
the infrared behaviour is consistent with a pole mass in 
the range 600 - 800 MeV. 

Our results give further support to the view that Nc ~ 
2 is a theory already close to the large Nc limit of QCD, a 
conclusion reached from the study of the string tension |^] 
and the glueball mass in SU{Nc)- 
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